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The structure of trapped binary mixtures of bosonic and fermionic atoms at zero temperature is 
studied using a modified Gross-Pitaevskii equation for the bosons which self-consistently includes the 
mean-field interaction generated by the fermionic cloud. The density of the fermionic component is 
described within the Thomas-Fermi approximation. The influence of the boson-boson and the boson- 
fermion s-wave interaction on the density profiles and the stability of the mixture is investigated 
systematically. Critical particle numbers for the mean-field collapse caused either by attractive 
boson-boson or by attractive boson-fermion interactions and for the onset of spatial component 
separation are discussed. It is shown that the interplay between the boson-boson and the boson- 
fermion interaction generates a rich and complex phase diagram. Finally, the specific properties and 
prospects of the boson-fermion mixtures available in present experiments are addressed. 
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I. INTRODUCTION 

Since the first realizations of Bose-Einstein condensa- 
tion in trapped dilute atomic gases [jl] this lively field of 
research has generated an amount of impressive exper- 
imental results illuminating basic quantum phenomena. 
Besides the studies using bosonic atoms first results were 
obtained on the cooling of fermionic atoms to a temper- 
ature regime where quantum effects dominate the prop- 
erties of the gas [|. One of the exciting prospects is 
the observation of a BCS transition of these degenerate 
Fermi gases to a superfluid state. 

Several successful attempts to trap and cool mix- 
tures of a bosonic and a fermionic species were reported. 
Quantum degeneracy was first reached with mixtures of 
bosonic 7 Li and fermionic 6 Li atoms ||, ||, [|. More re- 
cently, experiments to cool mixtures of different atomic 
elements, i.e. 23 Na and 6 Li @ as well as 87 Rb and 40 K 
H M, to ultra-low temperatures succeeded. These boson- 
fermion mixtures offer unique possibilities to study the 
effects of quantum statistics directly. Besides their inter- 
esting physics these systems can be used as an efficient 
tool to produce a degenerate single-component (spin- 
polarized) Fermi gas. Direct evaporative cooling is not 
applicable in a gas of spin-polarized fermions since s-wave 
interactions are suppressed by the Pauli principle. So the 
Bose gas, which can be cooled evaporatively, is used as 
a coolant. The fermionic cloud remains in thermal equi- 
librium with the cold Bose gas through boson-fermion 
interactions in the region where both clouds overlap. In 
this way the fermionic species is cooled sympathetically 
JTo| . Evidently, the shape of the two density distributions 
as well as their stability are extremely important for the 
success of the sympathetic cooling scheme. 

In this paper I want to investigate systematically the 
influence of inter-atomic interactions on the structure and 
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stability of binary boson-fermion mixtures at zero tem- 
perature. In extension to the work presented in the 
influence of boson-boson and boson-fermion interactions 
on the structure of the spatial density distributions and 
the stability of the system against collapse and phase 
separation is investigated. All combinations of attractive 
and repulsive interactions are considered and the impli- 
cations for present experiments are discussed. 

Previous studies of these mixtures mainly rely on the 
Thomas-Fermi approximation to treat the bosonic com- 
ponent Jl^, |l3) . This restricts the range of applicability 
significantly; the Thomas-Fermi approximation for the 
bosonic component is valid only for sufficiently strong 
repulsive boson-boson interactions and for large boson 
numbers. Other authors use a variational model, where 
the bosonic density profile is parameterized by a Gaus- 
sian |l4], [[5] . In many cases this does not provide enough 
flexibility to describe the density profiles accurately. 

To avoid these restrictions from the outset and obtain 
a model that can be applied for all values of the scat- 
tering lengths and for all relevant particle numbers the 
density of the bosonic component will be obtained from 
the self-consistent solution of a modified Gross-Pitaevskii 
equation which is coupled to the fermionic density dis- 
tribution through the boson-fermion interaction JTl| |l6| . 
The fermionic component is treated in Thomas- Fermi ap- 
proximation which actually is a very good approximation 
to describe a degenerate Fermi gas |l7], . 

The paper is structured as follows: In Section II the 
coupled set of differential equations that describes the 
bosonic and fermionic density distributions is derived and 
its numerical treatment is discussed. Using this universal 
tool the generic properties of binary mixtures are inves- 
tigated in detail in Sections III and IV. Depending on 
the particle numbers, the boson-boson scattering length, 
and the boson-fermion scattering length the instability of 
the mixture against mean-field collapse and spatial com- 
ponent separation is discussed. Finally, in Section |y| the 
consequences for the particular boson-fermion mixtures 
available in present experiments are summarized. 
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II. MODEL FOR DILUTE BOSON-FERMION 
MIXTURES 

We aim at the description of ground state properties of 
a dilute mixture of Ab bosonic and Af fermionic atoms in 
an external trapping potential at zero temperature. The 
atoms can be considered as inert interacting bosonic or 
fermionic particles, because the energy scales involved are 
sufficiently small, i.e., internal excitations of the atoms 
are not relevant. 



A. Mean-field approximation and effective 
Hamiltonian 

A simple starting point for the treatment of the quan- 
tum mechanical many-body problem of A — Ab + Af 
interacting bosons and fermions is the mean-field approx- 
imation. The many-body state that describes the binary 
mixture can be written as a direct product of a bosonic 
iVe-body state |$b) and a fermionic Ap-body state 
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|$b) ® |*f)- 
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Within the mean-field approximation each of these states 
is described by a symmetrized or antisymmetrized prod- 
uct of single-particle states. At zero temperature the 
bosonic state is a direct product of Ab identical single- 
particle states \4>b) 

(2) 



l$B> = \<h) 
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which is symmetric per se. The fermionic state neces- 
sarily is a Slater determinant, i.e. an antisymmetrized 
product of different single-particle states \ipi), \ipN F ), 



|* F > = A( |V>i) ® |^2) <8> ■ • ■ <8 \l/>N r )), 
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where A is the antisymmetrization operator. 

The mean-field picture is a reasonable starting point as 
long as correlations induced by the two-body interaction 
between the constituents are weak. The two-body po- 
tential that describes the interaction between the atoms 
consists of an attractive part at larger particle distances 
and a strong repulsion at short distances. Especially the 
short-range repulsion, the so called core, induces strong 
correlations between the atoms: The probability to find 
two particles within the range of the core is nearly zero, 
i.e., the two-body density shows a pronounced hole for 
particle distances smaller than the radius of the core. 
Correlations of this kind cannot be described in mean- 
field approximation. 

In order to use the mean-field approximation for the 
description of the interacting system one has to replace 
the full atom-atom potential by a suitable effective in- 
teraction that contains the relevant physical properties 
of the the original potential |l7| . Since the gases we are 
considering are extremely dilute and cold the particles do 
not resolve the detailed shape of the potential. Therefore 



one may replace the interatomic potential by a simple 
contact interaction for all partial waves. The strengths 
of the contact terms are adjusted such that the expecta- 
tion values of the contact interaction in two-body prod- 
uct states reproduce the two-body energy spectrum of 
the original atom-atom potential. This leads to a closed 
form of the effective contact interaction for all partial 
waves which involves only the scattering lengths of the 
original potential as discussed in detail in Ref. |l7| . 

Using the effective contact interaction one can set up 
the Hamiltonian of a binary boson-fermion mixture. It 
consist of a part Hb that involves only the bosonic com- 
ponent, a part Hp which describes only the fermionic 
component, and a part Hbf which represents the inter- 
action between the two species 



H = Hw + Hi + H 



BF- 



(4) 



The bosonic part contains the kinetic energy, the external 
trapping potential Ub(x) experienced by the bosons, and 
the s-wave interaction between bosons 
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where tub is the mass of the bosonic atoms and obo the 
s-wave scattering length of the boson-boson interaction 
(h = 1 throughout this paper). The formal distinction 
between bosonic and fermionic atoms is accomplished 
by projection operators Tlf and II? onto bosonic and 
fermionic single-particle states. Accordingly, the two- 
body projection operators n BB , IlfF, and II BF select 
pairs of two bosons, pairs of fermions, and boson-fermion 
pairs, respectively. 

The structure of the fermionic part of the Hamiltonian 
differs from that of the bosonic part: Due to the Pauli 
principle s-wave contact interactions between identical 
fermions do not contribute. The lowest order interaction 
contribution is therefore the p-wave contact interaction, 
which will be included here. The operator of the p-wave 
contact interaction is necessarily nonlocal, i.e., it involves 
the relative momentum operator = ~ (pi — j)j ) of the 
interacting particle pair. Thus the fermionic Hamiltonian 
reads 
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where ttif is the mass of the fermionic atoms, Up(x) the 
respective trapping potential, and czfi the p-wave scat- 
tering length |0 . 
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Finally, for the interaction between bosons and 
fermions, both, s-wave and p-wave terms are present 
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where tobf = 2m^mp/{mB + mp) is twice the reduced 
mass of a boson-fermion pair, obfo is the s-wave, and 
aBFi the p-wave scattering length of the boson-fermion 
interaction. 

In many cases the p-wave components of the interac- 
tions have only a marginal effect because the correspond- 
ing p- wave scattering lengths are small. Nevertheless, one 
can expect a significant influence of p-wave interactions 
on the structure and stability of the mixture provided the 
value of the p-wave scattering lengths is in the order of the 
s-wave scattering lengths or larger, e.g. in the vicinity of 
a p-wave Feshbach resonance |19j . The strong influence of 
p-wave interactions on the properties of purely fermionic 
systems was discussed in detail in Refs. ]l7], |l8| |. 

B. Energy functional 

Using the effective Hamiltonian ^ we study the 
ground state properties of boson-fermion mixtures at zero 
temperature. I adopt the language of density functional 
theory and construct, in a first step, an energy functional 
which connects the energy expectation value of the mix- 
ture with the one-body density distributions ub(x) and 
np(x) of the bosons and the fermions, respectively, 



-> E[n B ,n F }= d 3 x£(x) 
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For the construction of the energy functional the different 
parts of the Hamiltonian (Q) are considered separately 
and the energy density £(x) is decomposed into a purely 
bosonic part £b, a purely fermionic part £p, and a part 
£bf that contains the interaction between bosons and 
fermions 



£{x) = £b(x) + £f{x) + £bf(x)- 



(9) 



Using the many-body state (|J) the calculation of the 
expectation value of the bosonic part of the Hamiltonian 
Hb immediately leads to the well known Gross-Pitaevskii 
energy density pQj 
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The ground state density distribution n B [x) of the bosons 
is connected to the single-particle wave function (x\4>b) 



© by 
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which reflects the simple structure of the bosonic part of 
the many-body state. 

For the fermionic part the construction of the associ- 
ated energy density is not as straight forward. A general 
treatment of the fermionic many-body problem in mean- 
field approximation requires the solution of Hartree-Fock 
equations to obtain the single-particle wave functions 
p6[ . For the experimentally realized particle numbers of 
Np ?» 10 6 a self-consistent numerical treatment of these 
differential equations is not feasible. However, the large 
particle numbers allow the use of the Thomas-Fermi ap- 
proximation for the fermionic component fll?| . It is the 
lowest order of a semiclassical expansion of the energy 
density in terms of derivatives of the single-particle den- 
sity riF(x) All terms in the energy density that in- 
volve gradients or higher order derivatives of the density 
are neglected. To obtain the energy density in Thomas- 
Fermi approximation one calculates the energy expecta- 
tion value of the infinite homogeneous system using a 
Slater determinant of all momentum eigenstates up to the 
Fermi momentum kp with periodic boundary conditions. 
Next, one replaces the constant density np = fc F /(67r 2 ) 
therein by a density distribution np (x) . This leads to the 
fermionic contribution to the energy density |l7[ [l8| 
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where the first term represents the kinetic energy, the sec- 
ond term the external trapping potential, and the third 
term the p-wave fermion-fermion interaction. 

One can employ the Thomas- Fermi approximation also 
to calculate the contribution of the boson-fermion inter- 
action to the energy density 
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where the first term describes the s-wave and the second 
term the p-wave boson-fermion interaction. 

Notice that the Thomas-Fermi approximation affects 
only the contributions to the energy density associated 
with nonlocal terms of the Hamiltonian, that is the ki- 
netic energy and the p-wave interactions. For fermion 
numbers in the order Np » 1000 one can check explicitly 
that the Thomas-Fermi approximation is in good agree- 
ment with Hartree-Fock type calculations |lq| . Alterna- 
tively one can evaluate the next terms of the semiclassi- 
cal gradient expansion of the energy density 0. Both 
checks show that the Thomas-Fermi approximation gen- 
erally yields a very good description of the fermionic com- 
ponent for Np > 1000. 
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Some authors jlj, [L3| investigated the structure of 
boson-fermion mixtures using the Thomas- Fermi approx- 
imation also for the bosonic species. Here, however, one 
encounters the problem that the kinetic energy contri- 
bution of the bosons vanishes completely (not so for the 
kinetic energy of the fermions). This causes two stringent 
conditions for the applicability of the approximation: (a) 
the boson-boson interaction has to be repulsive, i.e., sys- 
tems with vanishing or attractive boson-boson interac- 
tion cannot be considered at all. (b) The contribution 
of the boson-boson interaction to the energy density has 
to be much larger than the neglected kinetic energy. For 
the investigation of the phase diagram as presented in 
the following sections this means a substantial and un- 
controlled limitation of the parameter range in which the 
approximation may be valid ]Tl| . Therefore, I employ 
the Thomas-Fermi approximation only for the fermionic 
component, where it is well suited for all types of interac- 
tions, and treat the bosons using the full Gross-Pitaevskii 
energy density. 



C. Euler-Lagrange equations 

The second step is the minimization of the total en- 
ergy E[n Bl n-p\ by functional variation of the boson and 
fermion density distribution. The variation shall be per- 
formed under the constraint of given numbers of bosons 
and fermions, N B and iVp. The constraints are imple- 
mented by two Lagrange multipliers, i.e. the chemical 
potentials /i B and fj,p, and an unrestricted variation is 
carried out for the Legendre transformed energy 



F[n Bl n F ] = £[ii B ,n F ] - ^i B N[n B ] - ^FN[n F ] 



(14) 



with a transformed energy density 

■F(cc) = £{x) — ^, B n B {x) — pl f tiy{x). (15) 

The values of the chemical potentials are chosen (a pos- 
teriori) such that the integrals over the boson (£ = B) 
and the fermion densities (£ = F) 

(16) 



N[n £ ] = l<Tx n^x) 



yield the desired particle numbers N B and Np, respec- 
tively. 

A necessary condition for the density distributions 
n B (x) — $ B (x) and n F (x) to describe a minimum of 
the transformed energy F[nB,n-p] is given by the Euler- 
Lagrange equations 
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One can insert the explicit form of J-(x), which involves 
the energy densities (|l0|), (|l2|), and (|l3|), and obtain two 



coupled equations that determine the ground state den- 
sities ns(f) and np(x). 

The Euler-Lagrange equation for the boson density re- 
duces to a modified Gross-Pitaevskii equation 
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In addition to the usual terms describing kinetic energy, 
trapping potential, chemical potential, and boson-boson 
interaction one obtains two contributions from the s- and 
p-wave interaction between bosons and fermions. The 
bosons thus experience additional mean-field potentials 
generated by the boson-fermion interactions that depend 
on different powers of the fermionic density distribution. 

For the fermions the Euler-Lagrange equation simpli- 
fies to a polynomial equation in n-p(x) because all terms 
in the energy density that involve the fermion density 
are local — a direct consequence of the Thomas-Fermi 
approximation 
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Again the boson density enters explicitly through the s- 
and p-wave boson-fermion interaction. The simultane- 
ous solution of Eqs. (|l8|) and ([19]) yields density profiles 
that describe a stationary point of the energy functional 
F[n B ,riY\- In principle one has to check explicitly that 
these coincide with a minimum of the energy. 



D. Numerical treatment 

On the basis of the coupled equations ( p"8| ) and ( 19j) 
the structure of binary boson-fermion mixtures is studied 
with special emphasis on the effects and the interplay 
between the boson-boson and the boson-fermion s-wave 
interaction. 

In order to reduce the number of physical parameters 
it is assumed that the p-wave scattering lengths are small 
and the corresponding p-wave terms are neglected. Thus 
the Euler-Lagrange equations for the boson and fermion 
density simplify to 
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for the boson density tib(x) = &%(x) and 
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for the fermion density. The index "0" which marked 
the s-wave scattering lengths is omitted here and in 
the following. Furthermore, we will restrict ourselves to 
parabolic trapping potentials with spherical symmetry 



TT-^X 
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where lj^ is the oscillator frequency and = (m^uj^)^ 1 / 2 
is the corresponding oscillator length for the bosonic 
(£ = B) and fermionic component (£ = F), respectively. 
In a magnetic trap the trap parameters of the fermionic 
species can be obtained from those for the bosons by 
scaling according to the different masses and magnetic 
moments. 

For the numerical treatment it is convenient to express 
all quantities with the dimension of length in units of the 
oscillator length Ib and energies in units of the oscillator 
frequency lob . Eventually, the physically relevant dimcn- 
sionless parameters to characterize the mixture are the 
particle numbers Nb and N-p, the ratios of the scatter- 
ing lengths and the oscillator length ob/^b and obf/^b, 
and the ratios of the masses tub / tof and of the oscillator 
lengths £b/£f- 

The solution of the coupled equations ( p0| ) and ( pi] ) 
is accomplished by an efficient imaginary time propaga- 
tion algorithm — also known as quantum diffusion al- 
gorithm p2j — which is embedded in a simple iterative 
scheme. A basic iteration step consists of two stages: (a) 
The fermionic density profile rip^x) is calculated from 
(2l|) with a chemical potential //p adjusted such that the 
desired fermion number Np is reproduced, (b) A sin- 
gle imaginary time step is performed using the mean- 
fields associated with the fermion density determined be- 
fore and the boson density of the previous iteration step. 
The resulting ub{x) is then normalized to the total bo- 
son number Nb and used as initialization for the next 
iteration step. In the beginning of the iteration cycle 
the boson density is initialized with a Gaussian profile 
chosen, e.g., according to a variational treatment of the 
isolated bosonic system. 

Typically 200 grid points are used to set up a radial 
lattice in the region where the bosonic and the fermionic 
density distributions overlap. The fermionic density pro- 
file outside the overlap region (where riB(f) = 0) is de- 
scribed using the analytic form (|2l]). The imaginary time 
propagation itself is implemented using cither a linearized 
time-evolution operator with a simple three-point dis- 
cretization of the radial part of the Laplacian or by a 
split operator technique employing a fast Fourier trans- 
formation. The size of a time step is chosen sufficiently 
small to guarantee the stability of the algorithm and the 



independence of the results. The convergence of the den- 
sity profiles is monitored during the evolution. 



III. ATTRACTIVE BOSON-FERMION 
INTERACTIONS 

In this section generic properties of binary boson- 
fcrmion mixture with attractive interactions between the 
two species (obf < 0) are discussed. Recent experiments 
indicate that mixtures of bosonic 87 Rb and fermionic 40 K 
atoms belong to this class of interactions ||, [|, [2^| . 

Throughout this and the following section it is assumed 
that the masses (m = tub — mp) and the oscillator 
lengths of the trapping potentials [l — Ib — £f) are 
identical for both components in order to reduce the pa- 
rameter manifold to a tractable size. 



A. Density profiles 

Some of the basic properties of mixtures with attrac- 
tive boson-fcrmion interactions become evident from the 
shape of the density profiles already. Figure |] shows the 
radial density profiles for mixtures with Nb = Np = 10 4 , 
cib/( — 0, and different strengths of the attractive boson- 
fcrmion interaction. Already the noninteracting case 
«bf/^ = (solid tlines) reveals a fundamental property 
of these mixtures: The fermionic density distribution has 
a much larger spatial extension than a bosonic distribu- 
tion with the same particle number. This is a direct con- 
sequence of the Pauli principle that forces the fermions 
to occupy "excited" single-particle states of the external 
potential which have a larger radial extension whereas 
the bosons in a Bose-Einstein condensate all occupy the 
ground state of the trap. This manifestation of the so- 
called Fermi pressure was observed experimentally § . 

The inclusion of an attractive boson-fermion inter- 
action leads to an enhancement of the boson and the 
fermion density within the central overlap region as it 
is shown by the broken curves in Fig. [l]. The bosonic 
profile narrows and the central density increases moder- 
ately. The effect on the fermionic component is more pro- 
nounced: within the overlap region the fermion density 
exhibits a high-density bump on top of the low-density 
profile. The central fermion density can easily be in- 
creased by a factor of four (dotted line in Fig. |l|). This 
structure can be understood by considering the mean- 
field potential experienced by the fermions: In addition 
to a shallow trapping potential the boson-fermion attrac- 
tion generates a tight potential well with the shape of the 
bosonic density distribution which causes a localized en- 
hancement of the fermion density profile. 

The boson-boson interaction has a strong influence on 
this density enhancement. The presence of a boson-boson 
repulsion broadens the bosonic distribution and reduces 
the maximum density significantly. This in turn strongly 
reduces the enhancement of the fermion density in the 
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FIG. 1: Boson (lower panel) and fermion density profiles 
(upper panel) for a mixture with Nb = Nf = fO 4 and 
o,b/£ ~ for different values of the boson-fermion scatter- 
ing length: cebf/^ = 0.0 (solid lines), —0.001 (dashed), and 
—0.002 (dotted). Both densities are given in units of the re- 
spective central densities ubo = Nb/(tv 3 ^ 2 £ 3 ) ~ 1796 t~ z and 
n F0 = (48iVF) 1/2 /(67r 2 ^ 3 ) « 11.7£~ 3 of the noninteracting 
system with the same particle numbers. 




FIG. 2: Boson and fermion density profiles for a mixture with 
N B = N F = 10 4 and a B /£ = 0.005 for different values of the 
boson-fermion scattering length: clbyII = 0.0 (solid lines), 
—0.01 (dashed), and —0.02 (dotted). The densities are given 
in units of the central boson and fermion density of the non- 
interacting system (compare Fig. |l|). 



overlap region. As is shown in Fig. g much stronger 
boson-fermion attractions are required to generate a sim- 
ilar increase of the fermion density in the presence of 
repulsive boson-boson interactions. However, the over- 
lap volume is much larger in these cases; thus a larger 
fraction of the fermions is contained in the high-density 
region. 



B. Simultaneous mean-field collapse 

For given particle numbers Nb and Np and given 
boson-boson scattering length ob there exists a maximum 



strength of the attractive boson-fermion interaction up to 
which the mixture is stable. Physically this stability limit 
is governed by the balance between the kinetic energy of 
bosons and fermions and the mutually attractive mean- 
field generated by the boson-fermion interaction. If the 
boson-fermion attraction becomes too strong — likewise 
if the boson or fermion numbers become too large — the 
attractive mean-field cannot be stabilized by the kinetic 
energy anymore and the mixture can lower its energy by 
increasing the boson and fermion density; both density 
distributions collapse simultaneously within the overlap 
region. 

In the numerical treatment of the coupled Gross- 
Pitaevskii problem the instability is indicated by a di- 
vergence of the central boson and fermion density during 
the imaginary time propagation. By monitoring the cen- 
tral density one can determine whether the systems is 
stable or collapses, i.e., whether the density converges or 
diverges for a given set of parameters JVb, Np, clb/£, and 
«bf/^- To determine the stability limit in presence of 
attractive boson-fermion interactions it is most useful to 
vary clbp/£ for fixed values of Nb, Np, and clb/1 until the 
onset of instability is reached. 

Figure || depicts the resulting stability limit in terms of 
a critical number of bosons N B as function of <zbf /i < 
for different ob/^ and Np. To obtain the curves I deter- 
mined the limiting value of asF ji that still allows a stable 
mixture for a set of boson numbers log 10 Nb = 1.5, 6. 
The maximum strength of the boson-fermion interaction 
of obf / ^ = — 0.05 shown in the plots corresponds to a 
scattering length of obf ~ — 1000aBohr ~ — 50 nm for a 
typical trap with £ — 1 /xm. 

The full curves in Fig. |^ show the critical boson num- 
ber for a mixture with vanishing boson-boson interac- 
tion (cib/1 = 0). Already a moderate boson-fermion at- 
traction with obf/^ = —0.005 causes a severe limitation 
of the boson number to Nb ^ 2300 in the presence of 
N-p = 10 5 fermions, as it is shown in panel (b) of Fig. |3[ 

The boson-boson interaction has a very strong influ- 
ence on the structure and stability of the mixture. If 
a weak boson-boson repulsion with a fraction of the 
strength of the boson-fermion attraction is included the 
system is stabilized significantly, i.e., the critical number 
of boson N B r is increased. This is depicted by the broken 
curves in Fig. |3| which correspond to different strengths 
as/^ > of the repulsive boson-boson interaction. For 
the previous example (czbf / ^ = —0.005, Np — 10 5 ) a 
boson-boson repulsion with ob/^ — 0.001 is sufficient to 
increase the maximum number of bosons in a stable sys- 
tem from N B r w 2300 for cib/^ = to far more than 10 6 . 
Moreover, repulsive boson-boson interactions can cause 
the absolute stabilization of the mixture against collapse 
induced by the boson-fermion attraction. That is, for 
given o,b/£ > and Np mixtures with |cibf/^| below a 
certain limiting strength are stable up to any number of 
bosons. This can be clearly seen in Fig. |(a): the bro- 
ken thick curves exhibit a kink and start to rise vertically 
at a specific value of cibf/£ which depends on ag/l and 
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-0.05 -0.04 -0.03 -0.02 -0.01 

a BF /£ 

FIG. 3: Critical number of bosons N§ above which the mix- 
ture collapses simultaneously as function of hbf/^ < 0. The 
thick curves correspond to repulsive boson-boson interactions: 
a B /e = (full line), 0.001 (dash-dotted), 0.005 (dashed), 0.01 
(dotted). Thin lines show attractive boson-boson interactions: 
as. 1 1 = -0.0005 (thin dash-dotted), -0.001 (thin dashed), 
—0.002 (thin dotted). The three panels represent different 
fermion numbers: (a) iV F = 10 4 , (b) N F = 10 5 , (c) iV F = 10 6 . 



Np. All systems with weaker boson- fermion attraction 
are stable against collapse for arbitrary iVe- 

An attractive boson-boson interaction, on the other 
hand, promotes the collapse of the mixture. It generates 
an additional attractive contribution to the mean-field 
experienced by the bosons which seeks to increase the 
boson density and has to be balanced by the kinetic en- 
ergy. Due to the boson-fermion attraction the fermionic 
component in turn feels a stronger mean-field attraction 
cause by the increased boson density. Therefore the crit- 
ical boson number is lowered by attractive boson- 
boson interactions as the thin curves in Fig. |^ show. 
One has to be aware that the character of the collapse 
changes if one reduces czbf/^ — * for fixed a^jl < 0. 



I 



For reasonably strong boson-fermion attraction the mu- 
tual mean-field couples the boson and fermion density 
strongly to each other such that both densities collapse si- 
multaneously in the overlap region. If the boson-fermion 
interaction is reduced this coupling weakens until bosons 
and fermions decouple for obf/^ = 0. In this case the 
boson density may still collapse due to the boson-boson 
attraction; the fermion density, however, is not affected 
and remains stable (see also Sec. IV C). 

The number of fermions in the system has a rather 
weak influence on the stability. The three panels in Fig. 
" correspond to (a) iV F = 10 4 , (b) N F = 10 5 , and (c) 
N-p = 10 6 . Mainly, the increasing number of fermions 
N-p enhances the effect of the boson-fermion attraction. 
Thus the stability limit for a given N-q is shifted towards 
smaller values of <2bf/^ if is increased. 



IV. REPULSIVE BOSON-FERMION 
INTERACTIONS 

In this section the generic properties of boson-fermion 
mixtures with repulsive interactions between the two 
species (a-QF > 0) are discussed. All 7 Li- 6 Li mixtures 
used in present experiments |Q, ||, |(| belong to this class 
of interactions. 



A. Density profiles 

Obviously a repulsive boson-fermion interaction will 
tend to reduce the overlap between the two species - 
in contrary to the attractive boson-fermion interactions 
discussed in the previous section. Hence the structural 
transition characteristic for this class of interactions is 
the spatial separation (or demixing) of the two species 
p"2| , |l6| ]. Figure ^ shows an example for the influ- 
ence of repulsive boson-fermion interactions of increas- 
ing strength on the density profiles of the mixture. The 
compact boson distribution in the trap center repels the 
fermionic cloud from the overlap region. With increasing 
boson-fermion scattering length obf /£ > the fermionic 
distribution gets more and more depleted until np(x = 0) 
reaches zero (see the dashed curve in Fig. fj). For even 
larger czbf/^ the fermions are expelled farther from the 
trap center and the overlap of both species is reduced 
continuously. The boson cloud in return is slightly com- 
pressed by the outer shell of fermions. 

The boson-boson interaction has a very strong influ- 
ence on the density profiles of bo sons a nd fermions. Sim- 
ilar to the case discussed in Sec. [II A a repulsive boson- 
boson interaction reduces the influence of the boson- 
fermion interaction. Much larger boson-fermion scatter- 
ing lengths are required to cause phase separation in the 
presence of repulsive boson-boson interactions. An ex- 
ample with as/i — 0.005 is shown in Fig. 0. Because 
the bosonic distribution is broadened due to the boson- 
boson repulsion, the fermion density is depleted in a sig- 
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cases there exists a second structural transition from the 
fermion-core configuration at smaller Nb to the usual 
boson-core configuration at large Nb ■ In deformed traps 
a rich variety of asymmetric configurations are possible, 



FIG. 4: Boson (lower panel) and fermion density profiles 
(upper panel) for a mixture with Nb = Nf = 10 4 and 
o,b/£ ~ for different values of the boson-fermion scatter- 
ing length: o.bf/1 = 0.0 (solid lines), 0.0007 (dash-dotted), 
0.0015 (dashed), and 0.003 (dotted). Both densities are given 
in units of the central densities of the corresponding nonin- 
teracting system as in Fig. |l|. 



0.08 



0.06 - 



: 0.04 



0.02 




FIG. 5: Boson (lower panel) and fermion density profiles 
(upper panel) for a mixture with Nb = A^f = 10 4 and 
a B /£ = 0.005 for different values of the boson-fermion scat- 
tering length: asF/f = 0.0 (solid lines), 0.005 (dash-dotted), 
0.015 (dashed), and 0.03 (dotted). The densities are given in 
units of the central densities ubo and wfo, resp., of the non- 
interacting system with the same particle numbers (compare 
to Fig. |l|). The dash-dotted curve for the boson density is on 
top of the solid curve. 



nificantly larger volume. A peculiar structure appears for 
interactions with clb/£ = clbf/( > 0: The fermion den- 
sity is nearly constant within the whole overlap region as 
the dash-dotted curve in Fig. | shows ||. 

Depending on the trap geometry the separated phase 
can exhibit different structures. If the fermionic species 
experiences a tighter confinement than the bosonic com- 
ponent (requires an appropriate combination of magnetic 
moments) an inversion of the separated configuration can 
appear: a central core of fermions surrounded by a thin 
boson shell which compresses the fermions. In these 



some of these were discussed in 16 



B. Spatial component separation 

In this section the influence of scattering lengths and 
particle numbers on the transition from overlapping to 
separated density distributions is quantified. Since the 
transition is smooth in most cases there is no unique def- 
inition for the onset of separation. As a gross criterion 
for separation I use the fermion density in the trap cen- 
ter and consider those configurations as separated that 
yield np(x = 0) = 0. Clearly, if the boson-fermion re- 
pulsion is increased then the overlap of the two density 
distributions decreases further. 

Figure |^ summarizes the findings on the character- 
istic number of bosons N^ p above which the fermion 
density vanishes in the trap center. The curves cor- 
respond to different values of the boson-boson scatter- 
ing length <2b /£■ For increasing boson-fermion repulsion 
cibf/£ > the characteristic boson number iVg° p drops 
rapidly. In the case a^jl — (solid lines) and Np = 10° 
a moderate boson-fermion repulsion with obf/^ = 0.005 
leads to phase separation if the number of bosons exceeds 
N^ p w 3300. 

The inclusion of repulsive boson-boson interactions 
leads to an interesting modification of the separation 
behavior. First of all the boson-boson repulsion stabi- 
lizes the mixture against phase separation, i.e., N^ p in- 
creases significantly with increasing a^/l > as the bro- 
ken curves in Fig. |6| indicate. Moreover, for sufficiently 
large fermion numbers the boson-boson repulsion causes 
a qualitative change of the dependence of N^ p on the 
boson-fermion scattering length. An example is shown 
in Fig. ^(c) for Np — 10 6 . Around a characteristic value 
of obf / £, which increases with increasing ob / 1, the limit- 
ing boson number for separation N^ p suddenly drops by 
one or more orders of magnitude and a steplike structure 
appears. For values of aBp/i smaller than the charac- 
teristic value associated with this step the boson-boson 
repulsion causes a substantial stabilization of the system 
against separation. Above the step the influence of the 
boson-boson repulsion on N^ p is weak. 

In addition to the particular behavior of N^ p the char- 
acter of the transition changes. One may regard the cen- 
tral fermion density np(0) as order parameter; finite val- 
ues of the order parameter characterize the overlapping 
phase, a vanishing order parameter indicates separation. 
In most cases, e.g. for as/^ = 0, the order parameter de- 
creases continuously with increasing Nb, i.e., the central 
fermion density decreases smoothly until the separated 
phase with np(0) = is reached. This behavior corre- 
sponds to a second order or cross-over transition. How- 
ever, for sufficiently large Np and ob/^, e.g., Np = 10 6 
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FIG. 6: Characteristic boson number iVg Cp for component 
separation as function of obf/^ > 0. The curves corre- 
spond to different values of the boson-boson scattering length: 
a B /e = (full line), 0.001 (dash-dotted), 0.002 (dashed), 
0.005 (dotted). The three panels represent different fermion 
numbers: (a) N F = 10 4 , (b) iV F = 10 5 , (c) 7V P = 10 6 . 



and as/i = 0.005 [see dotted curve in Fig. |(c)], the or- 
der of the phase transition changes. For values of obf/^ 
above the step the transition is of first order: With in- 
creasing Nb the order parameter decreases slightly until 
N^ p is reached, then it drops discontinuously to zero. 
For obf/^ below the value associated with the step the 
transition is still continuous which implies the existence 
of a tricritical point at the step. 



C. Collapse of the bosonic component 

If the repulsive boson-fermion interaction is supple- 
mented by a boson-boson attraction (as < 0) then a 
subtle competition between component separation and 
mean-field collapse occurs. 



For vanishing boson-fermion scattering length obf = 
the two species decouple. The bosonic component may 
undergo a mean- field collapse just like an isolated Bose- 
Einstein condensate if the boson-boson interaction is at- 
tractive. The critical particle number for this collapse of 
the bosonic component in a decoupled mixture can be pa- 
rameterized by Nft aBF=Q = 0.575 £/|cib|- This coincides 
with the result for the collapse a pure Bose gas po|, 24 
which was basically confirmed by experiment p5[ |26[ . 

One should notice that the character of the collapse 
induced by a boson-boson attraction differs from that of 
the collapse cause by an a ttrac tive boson-fermion inter- 
action discussed in Sec. [II B. In the latter case the 



boson-fermion interaction generates a mutual attractive 
mean-field that acts on both, fermions and bosons, and 
causes a simultaneous collapse of both density distribu- 
tions. In the case of a boson-boson attraction only the 
bosonic component experiences an attractive mean-field 
and may become unstable. For (Zbf > the fermionic 
density profile will be influenced by a collapse of the 
bosonic component but it will not collapse itself. The 
collapse of the bosonic component may, e.g., generate a 
collective excitation of the fermionic cloud. 

Although the collapse for a,B < concerns the bosonic 
density alone, the boson-fermion interaction modifies the 
stability properties significantly. We will concentrate on 
the case <zbf > 0; the collapse in presence of attractive 
boson-fermion interactions was discussed in Sec. Ill B 



already. The black curves in Fig. [?] show the critical 
number of bosons N B as function of the boson-boson 
scattering length as/^ for different values of the boson- 
fermion scattering length <zbf/^ > 0. The solid lines 
show the critical boson numbers for the decoupled mix- 
ture (obf = 0). 

The first general observation is that a boson-fermion 
repulsion destabilizes the bosonic component, i.e., the 
critical boson number is reduced if obf or N-p are 
increased. Physically this can be understood fro m the 
structure of the density profiles discussed in Sec. IV A. 



Due to the boson-fermion repulsion the outer fermionic 
cloud tends to compress the compact bosonic core. The 
boson density and thus the attractive mean-field are en- 
hanced which eventually promotes the collapse of the 
bosonic component. However, as can be seen from Fig. |^ 
the boson-fermion scattering length obf / ^ has to be sig- 
nificantly larger than |as/^| to have a noticeable effect. 

Evidently, component separation due to a boson- 
fermion repulsion has a large effect on the stability of the 
mixture with respect to collapse induced by the boson- 
boson attraction. The gray lines in Fig. |t] indicate the 
characteristic particle number N^ p for the onset of com- 
ponent separation for interactions with q,bf/^ = 0.01 
and 0.02. For weak attractive boson-boson interactions 
separation happens at lower particle numbers than the 
collapse of the bosonic component, that is N B r > N^ p . 
With increasing strength of the boson-boson attraction 
the critical particle number for collapse drops rapidly un- 
From that point on collapse and separa- 
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FIG. 7: Critical number of bosons N£ above which the 
bosonic component collapses as function of a-^/l < (black 
lines) as well as the characteristic boson number for com- 
ponent separation N^ p (gray lines). The curves correspond 
to different values of the boson- fermion scattering length: 
cibf/£ = (full line), 0.005 (dash-dotted), 0.01 (dashed), 0.02 
(dotted). The three panels represent different fermion num- 
bers: (a) iV F = 10 4 , (b) N F = 10 5 , (c) iV F = 10 6 . 



tion happen simultaneously, i.e. there is no stable sepa- 
rated configuration anymore. The corresponding curves 
for shown in Fig. exhibit a kink at this point and 
proceed towards larger |ae/^| like a smooth continuation 
of iVg Cp . There, as soon as component separation occurs, 
the bosonic cloud is compressed and the attractive mean- 
field due to the boson-boson interaction is amplified to 
such an extent that the collapse happens immediately. 



V. IMPLICATIONS FOR PRESENT 
EXPERIMENTS 

After investigating the generic properties of degenerate 
boson-fermion mixtures in the previous sections I now 
want to discuss the specific implications for the mixtures 
available in present experiments. In this section the mass 
difference and the different oscillator lengths of the trap- 
ping potentials for the two species are included explicitly. 



A. 7 Li- 6 Li mixture I 

Quantum degeneracy in a magnetically trapped di- 
lute boson-fermion mixture was first achieved by Tr- 
uscott et al. using a mixture of bosonic 7 Li atoms 
in the \F = 2,itif = 2) state and fcrmionic 6 Li with 
\F = 3/2, mp = 3/2). The same mixture was subse- 
quently used by Schreck et al. [|[ For this particular 
combination of states the boson-boson interaction is at- 
tractive and the boson-fermion interaction repulsive with 
scattering lengths |2j 



ob = — 1.46 nm , obf = 2.16 nm. 



(23) 



The attractive boson-boson interaction can cause the col- 
lapse of the bosonic component and leads t o a se vere lim- 
itation of the number of bosons (see Sec. HVCj ). This is 
a major hindrance for the implementation of an efficient 
scheme for the sympathetic cooling of the fermionic com- 
ponent. Ideally one would like to have a large bosonic 
cloud at a very low temperature to act as a coolant for 
the fermions. The attractive boson-boson interaction, 
however, restricts the number of particles in the Bose- 
Einstein condensate and hence sets a lower bound to the 
temperature of the bosonic cloud. In conclusion the low- 
est achievable temperature for the fermionic component 
is also restricted. 

The stability limits for this 7 Li- 6 Li mixture can be vi- 
sualized by a phase digram in the plane spanned by N-p 
and iVe as shown in Fig. ||. Each of the curves marks the 
limit of stability for a different size of the external trap- 
ping potential characterized by the oscillator length £-q 
for the bosonic species. The respective oscillator length 
£f for the fermions follows from simple scaling with re- 
spect to the different masses and magnetic moments. In 
the experiment of Truscott et al. the mean oscillator 
length is £b — 2.67 /mi [0, whereas Schreck et al. used a 
trap with £-q = 1.23 /jm H . For boson numbers above the 
stability limit the bosonic component is unstable against 
mean-field collapse. Narrowing of the trapping potential, 
i.e. reducing the oscillator length £b, enhances the effect 
of the interactions. 

Obviously the stability of the mixture depends mainly 
on the boson number JVb and the critical boson number 
is to a good approximation given by the corresponding 
value of the pure Bose gas. The fermionic component 
has essentially no influence on the stability because the 
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FIG. 8: Phase diagram for a 7 Li- 6 Li mixture with scattering 
lengths (p3|). The lines show the stability limit for different 
values of the oscillator length Ib of the trap as labeled. For 
Nb above the stability limit the bosonic component collapses. 



boson-fermion interaction is too weak. This is in agree- 
ment wi th th e experimental findings Q [3| . As discussed 
in Sec. IV C a stronger boson-fermion repulsion would 
destabilize the mixture anyway. 

Notice that the calculations presented here are re- 
stricted to isotropic traps. The prolate trap deforma- 
tion present in most experiments leads to a reduction of 
the critical particle number compared to a spherical trap 
with same mean oscillator length. For a pure Bose gas 
in a trap with uj z /uj r w 0.1 as in the experiment of Tr- 
uscott et al. B th e critical particle number is reduced by 
typically 20%~|28|. For the strongly deformed trap with 
uj z /Lu r w 0.02 used by Schreck et al. || this reduction is 
in the order of 40%. 



FIG. 9: Phase diagram for a 7 Li- 6 Li mixture with scattering 
lengths (||). The lines show the onset of component separa- 
tion for different values of the oscillator length Ib- For Nb 
above the limit the two species separate spatially. The gray 
curve shows the onset of separation for vanishing boson-boson 
interaction and Ib = \\jcm. 



the phase boundary at large N-p suddenly bend s tow ards 
smaller boson numbers. As discussed in Sec. IV B this 



structure is associated with a change in the order of the 
phase transition. At small Np or for shallow traps the 
transition is continuous, i.e. the central fermion density 
continuously goes to zero as Nb is increased. For large 
Np and tightly confining traps the transition is of first 
order, that is the central fermion density jumps to zero if 
the phase boundary is crossed. This effect is solely caused 
by the presence of the weak boson-boson repulsion. 



B. 7 Li- 6 Li mixture II 

The two Lithium isotopes can also be trapped in 
another combination of angular momentum states as 
was successfully demonstrated by Schreck et al. ||. 
The relevant s-wave scattering lengths for a mix- 
ture of 7 Li in the \F = l,mp = —1) state and 6 Li in 
\F = 1/2, tuf = -1/2) are both positive [|[ [To| 



a B 



0.27nm, aBF = 2.01nm. 



(24) 



Therefore this mixture is stable with respect to collapse 
but may undergo component sepa ration generated by the 
boson-fermion repulsion (see Sec. IV B ). 

The phase diagram of this mixture in the N^-Nb- 
plane is shown in Fig. ||. For values of Nb above the 
phase boundaries plotted for different oscillator lengths 
the mixture separates spatially. The onset of separation 
for the scattering lengths (|24|) depends in a nontrivial 
way on the particle numbers and the oscillator length. 
Especially for tightly confining traps with Ib < 0.15 fim 



We should like to point out that although the boson- 
boson scattering length is by one order of magnitude 
smaller than the boson-fermion scattering length it has 
a strong influence on the phase diagram also for shallow 
traps and small particle numbers. The gray curve in Fig. 
^ shows the phase boundary for £b — 1/im without the 
boson-boson interaction (compare with the black solid 
curve). The boson-boson repulsion stabilizes the system 
substantially, i.e. the onset of separation is shifted to 
much larger particle numbers. 

Schreck et al. H managed to produce such a mix- 
ture with Nb ~ 10 and Np « 2.5 x 10 4 in a trap with 
mean oscillator length of Ib = 1.0 /im (solid curve in 
Fig ||). According to our calculations (effects of the trap 
anisotropy are not included) the bosonic and fermionic 
density are still overlapping which is consistent with the 
experimental findings. However, an increase of the bo- 
son number by a factor 10 or a reduction of the fermion 
number to N-p « 10 3 would already lead into the regime 
of separation. 



12 




FIG. 10: Phase diagram for the 87 Rb- 40 K mixture with scat- 
tering lengths (p5|). The lines show the stability limit for 
different values of the oscillator length £b- For TVb above the 
limit the two species collapse simultaneously. The gray curve 
shows the stability limit for vanishing boson-boson interaction 
with £b = l^im. 



C. 87 Rb- 40 K mixture 

A first step towards an ultracold mixture of 87 Rb and 
the fermionic 40 K was reported by Goldwin et al. [9j and 
recently Roati et al. || managed to reach quantum de- 
generacy. This system is particularly interesting, because 
it shows a large negative boson-fermion scattering length 



ob = 5.25 nm , <jbf 



-13.8 nm. 



(25) 



Thus this mixture can undergo a simultaneous collapse 
of the bosonic and fermionic density distribution as dis- 
cussed in Sec. IHB . 

Figure ^ shows the corresponding phase diagram in 
the iVp-iVB-plane. One might expect that the strong 
boson-fermion attraction sets a very severe limit to the 
size of a stable mixture. In principle this expectation is 
correct as illustrated by the gray solid curve which cor- 
responds to the mere boson-fermion attraction for £b = 
1 /Ltm. Without the boson-boson interaction the mixture 
would collapse if the boson number exceeds « 300. 

Fortunately a moderate repulsive boson-boson inter- 
action is present which has a tremendous influence on 
the phase diagram. First of all the critical boson num- 
bers are increased by roughly four orders of magnitude. 
Moreover, below a certain fermion number marked by 
the vertical pieces of the curves in Fig. 10 the mixture is 
stable up to arbitrary boson numbers; this is the effect of 
absolute stabilization due to the boson-boson repulsion 
discussed in Sec. 



[II B 



Thus in a trap with £b — 1 M m 
any combination of boson and fermion numbers up to 
Nb ~ Nf ~ 10 7 would be stable although a strongly 
attractive boson-fermion interaction is present. 

The particular combination of scattering lengths of the 



87 Rb- 40 K mixture opens interesting perspectives for the 
implementation of a sympathetic cooling scheme. The 
strong boson-fermion attraction generates an enhance- 
ment of the boson and fermion density in the overlap 
region which is — due to the boson-boson repulsion — 
expanded over a large volume (see Fig. ||). Therefore 
a large fraction of the fermionic cloud overlaps with the 
bosonic distribution and enables a very efficient inter- 
species thermalization. At the same time the boson- 
boson repulsion stabilizes the system against simultane- 
ous collapse for all experimentally relevant particle num- 
bers. 



D. 



Na- Li mixture 



Recently, Hadzibabic et al. R reported the sympa- 
thetic cooling of a mixture of 23 Na and 6 Li to quantum 
degeneracy. To my knowledge there is no data on the 
inter-species scattering length for this system yet. How- 
ever, since a stable mixture with large numbers of bosons 
and fermions (N B w 2 x 10 6 , N F =» 1.5 x 10 5 ) was 
produced one can at least deduce bounds for the boson- 
fermion scattering length. Assuming ob = 2.75 nm |29 
one finds that the boson-fermion scattering length must 
obey <zbf > — 14.6nm in order to allow a stable mix- 
ture in this experiment. Moreover, the linear density 
profiles of the fermionic cloud shown in |?J suggest that 
the components are not separated. From the minimal re- 
quirement that the fermion density in the center is larger 
than zero one obtains an upper bound obf < 5.1 nm |p0| . 
Notice that these estimates do not account for the effects 
of the trap anisotropy and the finite temperature present 
in the experiment. 



VI. CONCLUSIONS 

In summary, I have presented a systematic study of the 
structure and stability of trapped binary boson-fermion 
mixtures at zero temperature. The bosonic component is 
described by a modified Gross-Pitaevskii equation which 
includes the mean-field interaction with the fermionic 
component. The fermionic density distribution is ob- 
tained within the Thomas-Fermi approximation. 

Depending on the strengths of the boson-boson and 
the boson-fermion interactions, characterized by the s- 
wave scattering lengths as and obf, respectively, differ- 
ent structural phenomena can be observed: 

as < : The collapse of the bosonic component domi- 
nates the phase diagram and sets a severe limitation 
to the maximum number of bosons in a stable mix- 
ture. The boson-fermion interaction has only small 
influence, however, attractive as well as repulsive 
boson-fermion interactions promote the collapse of 
the bosonic component. 



13 



ae > 0, obf < : The boson-fermion attraction en- 
hances the density of both species in the overlap 
region and can cause a simultaneous collapse of 
the bosonic and the fermionic cloud. A repulsive 
boson-boson interaction stabilizes the system, i.e. 
increases the critical particle numbers. 

«b > 0, cibf > : Repulsive boson-fermion interactions 
reduce the overlap of the two species and eventually 
lead to the spatial separation of the components. 
Again a repulsive boson-boson interaction stabilizes 
the mixture against separation. 

The rich variety of structures generated by the interplay 
of boson-boson and boson-fermion interactions has im- 
portant consequences for the experimental realization of 
sympathetic cooling schemes. As far as the density distri- 
butions are concerned efficient sympathetic cooling relies 
on (a) the mechanical stability of the mixture and (b) a 
large overlap of the two species. 

From this point of view mixtures with repulsive boson- 
boson interactions (<zb > 0) are superior since they do 
not suffer the severe limitation of the boson number im- 
posed by the collapse for as < 0. The boson-fermion 
cross section and thus the modulus of the boson-fermion 



scattering length should be sufficiently large to enable 
efficient inter-species thermalization. For obf > this 
may lead into problems with the separation of the two 
components, for <zbf < there may be a simultaneous 
collapse of both clouds. Fortunately, a moderate boson- 
boson repulsion can stabilize the mixture against both 
transitions. 

Regarding the structure of the density profiles a mix- 
ture with a B > and a B F < — like the 87 Rb- 40 K 
mixture discussed in Sec. V C — is especially appeal- 



ing for the implementation of an efficient sympathetic 
cooling scheme. It is stabilized against collapse by the 
boson-boson repulsion and exhibits a large overlap be- 
tween the species due to an extended boson distribution 
and a significantly increased fermion density within the 
overlap volume. 
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